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ABSTRACT 

Context. GJ 1214 is orbited by a transiting super-Earth-mass planet. It is a primary target for ongoing efforts to understand the emerg- 
ing population of super-Earth-mass planets around M dwarfs, some of which are detected within the liquid water (habitable) zone of 
their host stars. 

Aims. We present new precision astrometric measurements, a re-analysis of HARPS radial velocity measurements, and new medium- 
resolution infrared spectroscopy of GJ 1214. We combine these measurements with recent transit follow-up observations and new 
catalog photometry to provide a comprehensive update of the star-planet properties. 

Methods. The distance is obtained with 0.6% relative uncertainty using CAPScam astrometry. The new value increases the nominal 
distance to the star by ~10% and is significantly more precise than previous measurements. Improved radial velocity measurements 
have been obtained analyzing public HARPS spectra of GJ 1214 using the HARPS-TERRA software and are 25% more precise than 
the original ones. The Doppler measurements combined with recently published transit observations significantly refine the constraints 
on the orbital solution, especially on the planet's eccentricity. The analysis of the infrared spectrum and photometry confirm that the 
star is enriched in metals compared to the Sun. 

Results. Using all this new fundamental information, combined with empirical mass-luminosity relations for low mass stars, we 
derive updated values for the bulk properties of the star-planet system. We also use infrared absolute fluxes to estimate the stellar 
radius and to re-derive the star-planet properties. Both approaches provide very consistent values for the system. Our analysis shows 
that the updated expected value for the planet mean density is 1.6±0.6 g cirT 3 , and that a density comparable to the Earth (~ 5.5 g 
cm" 3 ) is now ruled out with very high confidence. 

Conclusions. This study illustrates how the fundamental properties of M dwarfs are of paramount importance in the proper charac- 
terization of the low mass planetary candidates orbiting them. Given that the distance is now known to better than 1%, interferometric 
measurements of the stellar radius, additional high precision Doppler observations, and/or or detection of the secondary transit (oc- 
cultation), are necessary to further improve the constraints on the GJ 1214 star-planet properties. 

Key words. Stars: individual :GJ 1214 — Astrometry — Techniques : radial velocities — Stars: late-type 



1. Introduction 

In the last few years, it has bec ome clear that a lar ge frac- 
tion of low mass stars (> 30% iBonfils et al.l |201 H) are or- 
bited by super-Earth-mass planets on relatively short period or- 
bits, including a handful detected in the l iquid water (habit- 
able) zones (e.g., GJ 581d a nd GJ 667Cc: iMavor et alJl2009t 
lAnglada-Escude et al.ll2012al respectively). These planets have 
small but non-negligible probabilities of transiting in front of 
their host stars. Given the favorable star-planet size ratio, we 
already have the technical means to begin spectroscopic char- 
acterization of their atmospheres. GJ 1214b was the first super- 
Earth found to transit in front of an M dwarf (Charbonnea u et al.l 
120091) and, as a consequence, it has received significant obser- 
vational and theoretical attention in the past three years. While 
the orbit is too close to the star to support hospitable oceans of 
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liquid water, GJ 1214b offers the first opportu nity to attempt at- 
mospheric characterization of a s uper-Earth dBean et al.l 12010: 
ICroll et alj201 UlBerta et aU 2012). but apparently contradictory 
results have emerged from these studies. Different gases in the 
atmosphere of a planet absorb light more efficiently at certain 
wavelength ranges. As a result, one should be able to measure a 
different e ffective transit de pth as a function of wavelength. For 
example. ICroll et all (1201 ll) reported excess absorption in the K 
band, which would be compatible with an extended and mostly 
transparent atmosphere with a strong absorb er in the near in- 
frared. On the other hand, iBean et all d2010h could not detect 
significant features in the optical transit depths. This result was 
confir med by the same group in a more ex tended wavelength 
range dBean etalJIioTlt l Desert et"aT]l201 ll) and has also been 
confirmed usi ng HST sp ectrophotometric observations in the 
near infrared (Bert a et al.ll2.Q12h . Such a flat transmission spec- 
trum favors a very opaque atmosphere with high concentrations 



1 



Guillem Anglada-Escude et al.: GJ 1214b revised 



of water vapor as the main source of opacity, indicating that the 
planet could be mainly composed of water. 

All these interpretations rely on atmospheric models that are 
strongly dependent on the planet's bulk properties, especially its 
mass, radius, mean density, surface gravity, and stellar irradia- 
tion. Prior to the detection of the planet candidate, GJ 1214 was 
a largely ignored M dwarf. As a result, some of its fundamen- 
tal properties had significant uncertainties (e.g., its distance). 
Uncertainties in the fundamental parameters of GJ 1214 prop- 
agate strongly into the planet's bulk properties, adding an extra 
element of uncertainty in discussions about the possible nature 
of its atmosphere. 

Probably the most significant measurement we provide here 
is the new measurement of its trigonometric parallax at 0.6% 
precisio n (previous parallax m easurement had a ~ 10% uncer- 
tainty, Ivan Altena et al.l [20011) . As a consequence, the lumi- 
nosity and mass of GJ 1214 have experienced significant up- 
dates as well. The orbital solution for GJ 1214 can also be up- 
dated using recently published transit observations and refined 
Doppler measurements obtained with our newly develop ed soft- 
ware (HARPS -TERRA, I Anglada-Escu de & But lerl2012l) . In ad- 
dition, the WISE catalog (ICutri et al.ll201 lb has also been re- 
cently released, adding four more absolute flux measurements 
of GJ 1214, thus enabling a comprehensive spectral energy dis- 
tribution adjustment and a more secure determination of T e ff. 
In Section [2] we present the new CAPSCam astrometric mea- 
surements and our re-analysis of public HARPS Doppler mea- 
surements. A new Keplerian solution for GJ 1214b is presented 
in Section 1231 In Section [3] we give an overview of the stellar 
properties in the light of the new distance measurement, its near 
infrared spectrum, and updated absolute magnitudes. Finally, 
Section [4] combines all of the transit observables with the new 
orbital solution and provides the posterior probability distribu- 
tions for the updated star-planet parameters. Our conclusions 
are summarized in Section[5] 

2. Observations and Data Reduction 

2.1. Astrometry and trigonometric parallax 

The trigonometric parallax of GJ 1214 has been obtained us- 
ing the Carnegi e Astrometric Planet Search Camera (CAPSCam, 
iBoss et al.l2009h installed in on the 2.5m duPont telescope of the 
Las Campanas Observatory (Chile). The observations span 26 
months (June 2009 to Sep 2012) and 10 epochs have been ob- 
tained at an average precision of 1 .0 milliarcsec (mas) per epoch. 
Each epoch consists of 20 or more exposures of 45 seconds 
each. GJ 1214 is significantly brighter than the average back- 
ground sources and would saturate the detector in less than 10 
seconds. However CAPSCam can re ad out a small p art of the 
array much faster than the full field (IBoss et al ] l2009h . For GJ 
1214, a window of 64x64 pixels is read out every 5 seconds, and 
all subrasters are added to the final full field image. The field 
of view of CAPSCam is 6.6 arcminutes wide and is typically 
rich in background stars, so a very robust reference frame with 
more than 30 objects can be used to correct for field distortions. 
Centroid extraction, source crossmatching, field distortion cor- 
rection, and astrometric solutions for all the stars in the field have 
been obtained using the ATPa astrometric software developed 
within the CAPS project (available upo n request). The m eth- 
ods and algorithms applied are outlined in lBoss et al.l d2009l) and 
Anglada-Escude et al. d2012bl) . 

A measurement of the parallax and proper motion requires 
fitting 5 astrometric parameters simultaneously: initial offset in 



Table 1. Local astrometric measurements of GJ 1214 
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R.A. and Dec, proper motion in R.A. and Dec, and the parallax 
itself. The formal uncertainties from a classic least-squares anal- 
ysis (e.g., from the diagonal of the covariance matrix) are usu- 
ally overoptimistic and do not properly account for correlations 
between parameters. To overcome these issues, we estimate the 
uncertainties in the astrometric parameters using a Monte Carlo 
approach. This is done by generating many realistic sets of mea- 
surements and measuring the standard deviation of the resulting 
derived parameters for the entire Monte Carlo-generated sample. 

To do this properly, a first realistic estimate of the epoch-to- 
epoch accuracy is needed. The outputs of the astrometric pro- 
cessing are the astrometric parameters of the target star and of 
all the other objects in the field. By combining the residuals of 
all the reference stars from all the epochs, we can compute the 
expected uncertainty per epoch. Doing this, we obtain an epoch- 
to-epoch accuracy of ~ 1.3 mas on both R.A. and Dec. Since the 
reference stars are fainter than the target, this is a conservative 
estimate of the real precision for the target. The star itself, which 
is not included in the reference frame, shows a standard devia- 
tion in the residuals of 1.0 mas/epoch, which is also consistent 
with the expected CAPSCam performa nce, assuming 2 minutes 
of on-sky observations per epoch (see lBoss et alJl2009h . Then, 
we simulate 10 5 synthetic sets of astrometric observations (same 
format as Table [TJ using the nominal parallax and proper mo- 
tion at the same epochs of observation. Random Gaussian noise 
with a single epoch uncertainty of 1.3 mas is injected into each 
synthetic data set and the 5 astrometric parameters are derived. 
The standard deviation of each parameter over the 10 5 solutions 
is the corresponding uncertainty. The obtained uncertainty for 
the differential parallax is 0.44 mas and corresponds to a rela- 
tive precision of ~ 0.6%. We note that this Monte Carlo method 
implicitly accounts for parameter correlation (see discussion in 
lFahertvetal1l2012l) . The differential astrometric measurements 
used to measure the differential parallax and proper motion of 
GJ 1214 are given in TableQ] The best fit solution over-plotted 
with the astrometric epochs is shown in Figure [1] 

The measured parallax is relative to the reference stars. 
Since they are not at infinite distances, these reference stars 
also have parallactic motion. As a result, the average paral- 
lax of the reference frame cannot be derived from the as- 
trometric observations alone. We obtain this average paral- 
lax (also called the parallax zero-point correction) using cat- 
alog photomet ry for the reference stars u sing the procedure 
described in Anglada- Escude et al] d2012bl) with 19 reference 
stars with good CAPSCam astrometry (RMS of the residuals 
below 2.0 mas/epoch) and reliable catalo g photometry (ff<18 
NOMAD and JHK S 2MASS photometn dZacharias et"aT]l2005b 
ISkrutskie et all2 006). For this set of observations, we determine 
a zero-point correction of -0.40 + 0.3 mas. This value, when 
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Fig. 1. Astrometric motion of GJ 1214 as a function of time. 
Top panels represent the motion in R.A.(left) and Dec. (right). 
Bottom panels show the residuals to the astrometric fit. The RMS 
of the residuals is below 1 mas. 

Table 2. Basic astrometric information and results from the anal- 
ysis of the astrometry. 



Parameter 




Value 


R.A. 




17 15 18.94' 


Dec. 




+04 57 49.7 1 


Catalog /4 A 


[mas yr 1 ] 


585 2 


Catalog /j Dec 


[mas yr~'] 


-752 2 


Relative yu^ A 
Relative ^ Dcc 




581.88 ±0.5 
-734.6 ± 0.8 



Relative parallax [mas] 69. 1 1 ± 0.4 

Zero-point correction [mas] -0.4 ± 0.3 

Absolute parallax [mas] 68.71 + 0.6 

Distance [pc] 14.55 ±0.13 

Notes. (1) 2MASS catalog dSkrutskie et al.ll2006h (2) LSPM-NORTH 
catalog lLepine & Sharal d2005f) 

added to the relative parallax, produces the final absolute paral- 
lax measurement in Table [2] The inverse of the parallax in arc- 
seconds is the distance in parsecs (pc) and amounts to 14.55 ± 
0.13 pc. The updated BVRIJHK S , Wl, W2, W3 and WA absolute 
magnitudes for GJ 1214 are presented and discussed in Section 

m 

We note that the measured proper motions are also differ- 
ential and contain an unknown offset due to the unknown aver- 
age motion of the background stars (e.g., all galactic plane stars 
move roughly in the same direction). Even though catalog val- 
ues are less precise than the ones we obtain, catalog proper mo- 
tions are usually corrected for proper motion zero-point ambi- 
guities. Differential measurements as well as the suggested val- 
ues for the prop er motion of GJ 1214 (LSPM-NORTH catalog, 
iLepine & Sharal2005l) are also provided in Table [2] 

2.2. Radial velocity measurements 

The ESO public archiv qj contains 21 reduced HAR PS spectra 
of GJ 1214 obtained by ICharbonneau et all d2009l) . We rean- 
alyzed these public HARPS spectra using our software called 
HARPS -TERRA (HARPS T emplate Enhanced Radial veloc - 
ity Re-analysis Application, Anglada-Escude & Butlerl l2012h . 
HARPS-TERRA derives differential RV measurement by con- 
structing a high signal-to-noise ratio template from the obser- 
vations and matching it to each spectrum using a least-squares 
approach. GJ 1214 is faint at optical wavelengths (V~14). As a 
consequence, very long exposures (2400 sec) were required to 
obtain any signal at all. Still, the typical S/N at 6100 A is only 
~10 so careful construction of the template is a key element in 

1 http : / /archive . eso . org/wdb/wdb/eso/repro/ form 



Table 3. HARPS-TERRA radial velocity measurements. 
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order to achieve the maximum precision. The standard setup of 
HARPS-TERRA for M dwarfs uses all the echelle apertures red- 
der than the 22nd one (4400A< A < 6800A) and adjusts a cubic 
polynomial to correct for the variability of the blaze function 
across each echelle order. Detailed algorithms and performance 
of HARPS-TERRA on a rep r esent ative sample of stars is given 
in Anglada-Escude & B utlerl (120121) . It is wo rth noticing that er- 
ror ba rs listed in [3] are smaller than those in Charbonneau et alj 
( 2009), suggesting a more optimal usage of the Doppler informa- 
tion in the spectra. However, reduced error bars do not guarantee 
a more accurate orbital fit when the RV variability is dominated 
by systematic noise (either instrumental or stellar). To account 
for this, we included an unknown noise term in our Monte Carlo 
Markov Chain runs which had a substantial effect on the derived 
distributions, especially on the orbital eccentricity (see discus- 
sion in Section l2~3l ). 

2.3. Orbital solution 

As discussed in ICarter et"aTI d201 ll) . uncertainties in the orbital 
parameters (especially the eccentricity) are the main limitation 
in the derivation of precise star-planet parameters. Therefore, 
we are interested not only in the favoured orbital solution, but 
also in obtaining a realistic numerical representation of the pos- 
terior density function of the parameters to propagate them 
into the estimates of the star-planet properties. These samples 
are generated using a Bayesian Monte Carlo Markov chain 
(MCMC) method. We use custom-made MCMC software to 
combine RV with transit observations and to obtain such dis- 
tributions. The softwa re is based on the MCMC tools used in 
lAnglada-Escude et alj d2012b ) adapted to the Doppler plus tran- 
sit problem, and uses a Gibbs sampler. Jump scales of the Gibbs 
sampler are initialized using the diagonal elements of the co- 
variance matrix at the maximum likelihood solution, and they 
are adjusted in the first 10 s steps (the "burn-in" period) to accept 
between 10% and 30% of the jump proposals for each param- 
eter. As discussed below, we assume flat priors for all the sam- 
pling parameters. As a consequence, the probability of accepting 
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a new state is just the ratio of likelihoods at the current and pro- 
posed jump position. An outline of the general MCMC meth od 
applied to the Keplerian problem can be found in lFordl d2006h . 

The Keplerian model for the ra dial velocities is the same as 
in Anglada-Es cude & Tuomil d2012l) and the transit observations 
are predicted using the recipes given in the NASA Exoplanet 
Archive webpage 0. Transit observations strongly constrain the 
orbital period, but also put strong restrictions on the sum of 
the initial mean anomaly Mo, the argument of the node to (also 
called, mean longitude Aq = Mq+oS), and the product e cos to (ec- 
centricity and argument of the node). Accordingly, we use Aq as 
a free parameter. Since the transit instants explicitly depend on 
e cos to, e sin to and e cos to are sometimes used as the MCMC 
sampling parameters. We also tested this approach, and, while 
it has some desirable properties (e.g., symmetry around zero 
eccentricity)- this parameterization imposes an implicit priori 
that sever ely biases the estimates of e to artificially larger values 
(e.g., see lFordl2006l;lBarros et al.l201ll for similar discussions). 
Therefore, at the sampling level, we choose those free parame- 
ters that, from our point of view, should have flat prior distribu- 
tions. These are: orbital frequency 1/P, RV semi-amplitude K, 
mean longitude Aq, argument of the node to, and eccentricity e. 
Despite to being poorly defined and uncertain at low eccentrici- 
ties, the convergence of the MCMC was substantially improved 
thanks to the use of Aq instead of Mo, so we find this parametriza- 
tion convenient and sufficient for our purposes. All the approxi- 
mate parameter values are already known from previous studies 
and, therefore, we initialize the MCMCs wi thin 3 standard de- 
viation s of the orbital solution proposed by ICharbonneau et all 
(120091) . 

Real uncertainties in radial velocity measurements are diffi- 
cult to estimate properly, especially when the star is active. To 
account for that, we model the uncertainty of the z'-th measure- 
ment as cr? = e? + s 1 , where e,- are the formal uncertainties de- 
rived from HARPS-TERRA (third column in Table [3}, and s is 
a jitter parameter that accounts for the extra systematic noise. 
In this context, the jitter parameter s is treated as any other free 
parameter. 

We assign a constant uncertainty to each transit instant of 
400 sec. This uncertainty is 10 times larger than the typical re- 
ported formal errors but it is more consistent with the differ- 
ent instants of transit measurements obtained simultaneously by 
different groups (see Table [8). While the MCMC convergence 
properties are greatly improved, 400 sec is a small fraction of 
the time-baseline covered by the transits and, as a consequence, 
the period is still strongly constr ained. The t r ansit t i me measure- 
ments have been extracted fr om Sada et all (1201 Oh . ICarter et aT] 
(1201 lh . and lKundurthv et all (l201lh . 

In table |4] we provide the expected values for the most rel- 
evant combination of parameters for physical and transit pre- 
diction computations as derived from the combination of 10 
MCMC runs with 10 7 steps each. Figure [2] is an illustration 
of the maximum a posteriori probability (MAP) solution and 
is given only for illustration purposes. Since the amplitude of 



2 http://exoplanetarchive.ipac.caltech.edu/applications/TransitSearch/ 
guide/algori thms.html 

3 The implicit prior comes from the fact that x = e sin to and y = 
e cos to are cartesian coordinates derived from the polar radial coordi- 
nate e and the angle to. Because MCMCs are based on the integration of 
a probability distribution, to preserve the flat prior choices for e and to, 
one would need to include the Jacobian of the transformation (1/e) as 
a prior when using x and y as the Markov chain sampling parameters. 
The use of such a prior is undesirable due to numerical stability issues 
(divergence for e = 0). 




-0.4 -0.2 0.2 0.4 
Orbital phase [days] 

Fig. 2. HARPS-TERRA radial velocity measurements folded to 
the orbital period. The instant of transit is depicted as a dashed 
vertical line. 
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Fig. 3. Left. Distribution of MCMC states for the derived quan- 
tities ecosw and e sin to as obtained from 1 MC chain of 10 7 
steps (red) and the combination of 10 chains of 10 7 steps each 
(black). Only 1 % of the steps are included in this plot to improve 
visualization. Right. Marginalized probability distribution of the 
eccentricity is shown in black (arbitrary normalization). The cor- 
responding cumulative distribution (red histogram) shows that 
eccentricities higher than 0.23 are ruled out at a 95% confidence 
level. 



the signal is small compared to the noise, the eccentricity is 
still poorly constrained. The posterior distribution of the eccen- 
tricity has a maximum very close to and monotonically de- 
creases with e (see Figure |3]l implying that close to circular 
orbits are favored. Because of their significance for the transit 
observables (e.g., see Section 0), an illustration of the MCMC 
samples of the derived parameters e cos to and e sin to is also 
provided in figure [3] Previous estimates of the e c centri city of 
GJ 1214b ICharbonneau et all (120091) : ICarter et al.1 (1201 ll e.g.. ) 
seemed to favor eccentricities close to 0.1, but still compatible 
with 0. In initial tests, we recovered this same slightly higher 
eccentricity value if the noise parameter s was fixed to 0. This 
indicates that previous estimates of the orbital elements were bi- 
ased due to a non-consistent treatment of the RV uncertainties. 
Even though stellar jitter dominates the error budget {s ~ 2.9 
m s _1 ), the RMS of the MAP solution using the new HARPS- 
TER RA measurements is low er (3.4 m s _1 ) than the one reported 
bv ICharbonneau et al.l (|2009, 4.4 m s _1 ). Given the strong ef- 
fect of the jitter parameter s on the result, and until more RV 
measurements become public, we strongly recommend the use 
of the updated solution in Table [4] for any future work (e.g., in 
searching for seconda ry transits). Applying the s ame 95% con- 
fidence level used by ICharbonneau et all ^2009), we obtain an 
upper limit to the eccentricity of 0.23. 
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Table 4. Expected values of the most useful parameter combi- 
nations obtained using RV and transit observations. Numbers in 
parentheses represent the standard deviation of the distributions 
(last two significant digits of the expected values). 



Parameter 



P [days] 


1.580400(14) 


Kims" 1 ] 


10.9(1.6) 


<k [deg] 


210.9(6.2)" 


ecos to 


-0.033(55) 


esinai 


-0.044(90) 


Jitter s [m s -1 ] 


3.6(1.1) 


y 


0.35(1.2) 


Other parameter combinations 




e 


<0.23" 


Msini [Mju P ] 


0.0195(28)" 


M sin i [M s ] 


6.20(91) c 


a[AU] 


0.0148 r 




^transits 


30 




21 


RMS RV [ms-»] 


3.5 





Notes. {a) Julian date of the reference epoch on which A is computed 
is the first epoch of the RV data at T Q = 2455036.57372 days. For 
a circular orbit A is equal to the mean anomaly M at the reference 
epoch. <A) Distribution of the eccentricity peaks close to (see Figure 
FJ). w Assumes M»=0.176 M . The uncertainty in the mass of GJ 1214 
is considered in the final proposed parameters for the star-planet system 
given in Section|4] 



Table 5. Updated absolute photometr y for GJ 1214. Johnsons - 
Cousins BVRI photometry is from dDaws on & Forbes 1992). 
JHKs photometry is from the 2MASS catalog (Skrutskie et al. 
2006). The Wl, W2, W3 and W4 are the mid-infrared bands 
from the WISE Preliminary Data Release dCutri et all 1201 ll) 
(central wavelengths are 3.4 yU, 4.6/im, 12/mi and 22/im). 
Properties derived from the fits of the absolute photometry to 
the BT-Settl-2010 model grid are also provided. 





LJ.JJZ, ± U.UoU 


M,- 
Ivly 




LVLR 


1 ? 44 1 +0 04 ^ 


M; 


11.608 ± 0.042 


Mj 


8.934 ± 0.041 


M„ 


8.274 ± 0.041 


M Ks 


7.964 ± 0.038 


M m 


7.781 ±0.041 


M w2 


7.614 ± 0.039 


Mm/3 


7.407 ± 0.041 




7.204 ±0.178 


Teff [K] 


3252 ± 20 1 


L, [10- 3 L G ] 


4.05 ±0.19 


R, [«q] 


0.201 ± 0.010 


[Fe/H] ,,/,„, 


+0.13, +0.05 2 


[Fe/H] spec 


+0.20 3 


\M/H\ SP ec 


+0.15 3 



Notes. (1) Uncertainty in T e ff o nl y acc oun ts for statis ti cal er - 
ro rs. (2) | Schlaufman & Lauehlml hoiOh : iNeves et alj d2012h 
(3) lRoias-Avala et al.1 d2012h 



Another interesting parameter is the minimum mass of GJ 
1214b. Even though the new RV amplitude K is smaller (~ 
10.9 m s" 1 ) than the previously reported one (~ 12.5 m s~'), 
the updat ed minimum mass d oes n ot change significantly com- 
pared to Char bonneau et all d2009l) due to the similar relative 
increase in the updated stellar mass derived from the new dis- 
tance determination (see Section 0J. The origin of the slightly 
smaller amplitude is likely caused by stellar activity. By analyz- 
ing D oppler measuremen ts on known active M dwarfs (e.g., AD 
Leo, iReiners et all 1201 2l) . we found that the classic HARPS- 
CCF approach always produces a larger Doppler amplitude than 
the one derived using HARPS-TERRA (between 1 . 1 to 1 .5 times 
larger) if a candidate signal is caused by activity. The likely ex- 
planation for this is the different sensitivities of the two meth- 
ods to the changes in the stellar line-profiles. While one could 
use this effect to obtain a further diagnostic to check the reality 
of low-mass candidate planets (under investigation), in this case 
it means that the most likely source of systematic noise comes 
from the star rather than from the instrument. Doppler follow-up 
of the star is needed to better understand the origin of the extra 
noise and to further constrain the orbital solution. 

3. Updated properties for GJ 1214 

3.1. Metallicity 

Two independent tech n iques, the photometric method by 
Schlaufm an & Laughlinl (l2010h and the K -band spectroscopic 
[Fe/H] index by iRoias-Avala et alJ (I2012I) . agree on the metal- 
richness of GJ 1214, with [Fe/H] = +0.28 dex and [Fe/H] = 
+0.20 dex, respectively. However, since the [Fe/H] photomet- 



ric calibrations depend on the distance of the star (to obtain 
M^), the photometric value of GJ 1214 needs to be re-calculated 
using its updated distance. The updated absolute magnitudes 
are given in Table [5] The updated photometric [Fe/H] values 
are [Fe/H] = +0.13 dex and [Fe/H] = +0.05 dex, using the 
ISchlau fman & Laug hlinl (l2010h calibration and the recent cali- 
bration bv lNeves et al.l (|2012), respectively. Considering the dis- 
persions associated with the calibrations (cr 0.1-0.15 dex), the 
estimates are consistent with each other, corroborating that GJ 
1214 is a solar or super-solar [Fe/H] star. 

A comparative approach can also be performed to confirm 
the metal-richness of GJ 1214. Figure|4]shows the K-band spec- 
tra of Gl 699 (Barnard's star), Gl 231. IB, and GJ121/Q. Gl 699 is 
the second nearest M star to the solar system and its kinematics, 
Hq, activity, and [Fe/H] meas urements are consistent with Gl 69 9 
being an old disk/halo star (iGizisI 1 997 : Roias-Av ala et al.l2012l) . 
Gl 231. IB is the low-mass compa nion of a nearly solar me tallic - 
ity GOV star ([Fe/H]=-0.04 dex. lValenti & Fischerl d2005h '). The 
overall shapes of the K-band spectra of all three stars are quite 
similar (same spectral type), and the morphology of the surface- 
gravity-sensitive CO bands corresponds to M dwarf stars (log 
g ~ 5). However, all the absorption features of GJ 1214 are 
stronger tha n the ones exhibited by m etal-poor Gl 699 ([Fe/H] 
= -0.39 dex. lRoias-Avala et all d20 121) 1 and solar metallicity Gl 
23 LIB. Therefore, also in a relative sense, the strong absorption 
features favor a high metal content for GJ 1214. 



4 These spectra are part of the K-band spectral atlas by 
IRoias-Avala et alj i2012h . and available to the community in the online 
version of that article. 
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Fig. 4. K-band spectra of GJ 1214 (top), Gl 231. IB (middle), 
and Gl 699 (bottom). All the stars have the same K-band spectral 
type (same overall shape of their spectra), but different metallic - 
ities, as indicated by the strengths of their absorption features. 
GJ 1214's metallicity should be at least equal to or higher than 
solar. 



3.2. Stellar temperature, luminosity and radius from models 

We first provide an overview of previous methods and de- 
terminations of the temperature, luminosit y, and radius of GJ 
1214 . In the discovery paper of GJ1214b, (Charbonneau et alj 
20091) estimated a T eff = 3026K and a R*= 0.21 IRq for GJ 1214. 
Kundurthv et "all (1201 ll) obtained T eff =2949K and R* = 0.21 1R Q . 
Constraining param e ters to stell a r evo lution isochrones by 
iBaraffe et all dl998l) . ICarter et all d2011l) led to estimates of 
T eff =3170K and R*= 0.179R o , assuming that GJ 1214 is an "old 
star". However, all thes e results were based on th e distance deter- 
mination of ~ 13 pc by Ivan Altena et all d2001l) . and reveal that 
the effective temperature of GJ 1214 is quite uncertain, given 
all the degeneracies in the evolutionary model used and a low 
precision distance estimate. 

We attempted to get an estimate of the stellar radius using 
up-to-date interfero metric empirical calibrations. For example, 
iKervella et al.l d2004l) provides an empirical luminosity-radius 
calibration that extends to low mass stars and uses absolute mag- 
nitudes as the only inputs. We tried different combinations of 
photometric bands obtaining inconsistent results. All the rela- 
tions involving optical bands (BVRI) gave values larger than 0.4 
Rq, which are very unrealistic considering the spectral type of 
GJ 1214. The relations restricted to infrared colors(7, H, and K) 
provided estimates a bit more realistic (between 0.15 to 0.2 R ) 
but still not fully comp atible with each other . We also tried the 
calibration provided by De morv et alj d2009l) that used K band 
photometr y to minimize the m etallicity effect on optical magni- 
tudes. In iDemorv et al.l d2009l) . new measurements of 6 new M 
dwarfs were presented and a new radius-luminosity calibration 
was derived. The y found a remarkable agreement with the evolu- 
tionary models of Baraffe et al. ( 1998) if the measured radii was 
plotted against the K absolute magnitude. Using this approach, 
we obtain a stellar radius of 0.193 R which, at least, seems to 
be in the expected range. 

Given that the state-of-the-art empirical relations are not en- 
tirely self consistent, we also obtained a new estimate of the ef- 
fective temperature, luminosity, and radius of GJ 1214 using up- 
dated BVRIJHK S W1 W2 absolute photometry (Tabl ed together 
with the model atmosphere grid BT-Settl-2010 dAllard et alJ 



1201 llPl The lack of indicators of youth in GJ 1214, such as H a 
emission, along with i ts space motion, giv es an age estimate of 
3-10 Gyr for the star dReid & Gizisll2005l) . This age estimate is 
also s upported by its long rotation period (Charbonneau et al. 
2009). Therefore, we fixed the surface gravity of the syn- 
thetic spectral models to log g = 5.0 and [M/H]=0.0. Figure [5] 
shows the bolometric luminosity as function of absolute mag- 
nitude for the BT-Settl-2010 grid. The inferred bolometric lu- 
minosities with the JHKsWlW2 photometry are consistent with 
each other, and higher than the luminosities obtained with the 
BVR I photometry. Pr evious synthetic models (e.g. NextGen; 
lHauschildt et al.lfl999l) showed a lack of flux in the K band when 
compared with observed spectra of low-mass stars. New solar 
abundances and the inclusion of dust grain formation seems to 
have solved most of the previous discrepancy, allowing the BT- 
Settl-2010 models to reproduced fairly well the infrared sp ectral 
energy distribution (SED) of M-dwarfs dAllard et al.ll2011l how- 
ever, the FeH opacity data is still incomplete for this region). 
Although models have also been improved at shorter wave- 
lengths (BVRI), they remain too bright in the ultraviolet and vis- 
ible part of the M dwarf spectra, pos sibly due to missin g sources 
of opacity in the modeling process dAllard et alJl201 ll) . The ef- 
fective temperatures and radii corresponding to the best SEDs 
fit to the 9 wavebands BVRIJHK S W\ W2, only BVRI, and only 
JHK S W1W2, obtained from their respective bolometric lumi- 
nosities in Figure [5] are shown in Figure [6] Given that 1) bolo- 
metric luminosities obtained with JHK S WIW2 photometry are 
consistent with each other, 2) the models provide a better fit at 
longer wavelengths, and 3) the empirical luminosity-radius pro- 
vides more self-consistency using nIR colors, the fit with only 
JHK S W1W2 provides the most reliable results and is the one to 
be used in deriving further properties of the star-planet system. 

The SED fitting using only these nIR bands is further sup- 
ported by other results. For example in Section |H we will ob- 
tain an independent stellar radius estimate (0.211 ± 0.011 R ) 
combining the stellar mass with direct observables from the 
light curve and Doppler data. Such an estimate is compatible 
with the value we discussed in this section (~0.2 R ). Also, 
the estimate of the effective temperature (T e ff ~ 3250/f) is 
in excellent agreement with the effe ctive temperature derived 
from water absorption in K-band by Roias-Av ala et al.l d2012i 
T e ff =3245 K). Furthermore, next Section shows this T e ff also 
agrees w ith the one derived from SED fits to the evolutionary 
models dBaraffe et al.lll998l) . 

3.3. Stellar mass 

Although most of the previously reported mass-radius estimates 
of GJ 1214 agree, we not e that they all assume th e trigonomet- 
ric distance estimate from lvan Altena et al.l (12001). This distance 
determination has an uncertainty of 10% which is not always ac- 
counted for in the aforementioned predictions. We also note that 
the Ml V-M4V spectral type range spans a broad range of stellar 
masses (from 0.4 to 0.15 M ) and effective temperatures (3500 
to 2800 K). This makes any mass estimate very sensitive to un- 
certainties in the parall ax mesurement . Given the updated dis- 
tance, we can apply the Del fosse et al.l (l2000l) (hereafter DF00) 
relations to derive a new mass for GJ 1214 using the absolute J, 
H and K magnitudes. The masses derived from each photomet- 
ric band are M (y) = 0. 1 74 M ,M (H) = 0. 1 77 M and M w =0.177 
o, which are in very good agreement. 



5 http://phoenix.ens-lyon.fr/simulator/index.faces 
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Fig. 5. Luminosity as function of absolute magnitude for all the 
wavebands listed in Table [5] Open diamonds represent the pho- 
tometry of GJ 1214. The color asterisks represent the BT-Settl- 
2010 values, and the color dotted lines the linear interpolations 
between them. The dashed line indicates the mean bolomet- 
ric luminosity estimated using all the photometry (log(L/L )= - 
2.57). The dotted-dashed line indicates the mean luminosity with 
only the BVRI photometry (log(L/L )= -2.73), and the 3 dotted- 
dashed line the adopted mean luminosity for GJ 1214, derived 
using only the JHK S W\W\ absolute magnitudes (log(L/L )= 
-2.44). The lowest bolometric luminosity corresponds to the I 
magnitude, where the SED is very steep. 



Fig. 6. Spectral energy distribution fits for GJ 1214 for the ef- 
fective temperatures and radii derived from photometry. The 
magenta, green, and blue SEDs are the spectral templates that 
represent the values obtained with only BVRI photometry, only 
JHK S WIW2 photometry, and the 9 wavebands, respectively. All 
spectral templates have solar metallicity, log g = 5.0 and distance 
equal to 14.47 pc. The photometry data of GJ 1214 are plotted 
as black dots. Color dots indicate flux levels of each spectral 
template integrated over the corresponding filter bandwidth, de- 
picted by the horizontal black lines. The spectral template with 
T C //=3250K and R=0.2R o provides the best fit to the infrared 
photometry of GJ 1214. 



Another way to obtain an estimate of the mass for GJ 
1214 is to compare the absolute magnitud e with the s y ntheti c 
colors from th e evolutionar y mode ls in iBaraffe et al.l (Il998h . 
Unfortunately, Baraffe et al. (1991) does not contain models 
for stars with [Fe/H]>0.0. As a result the optical colors (e.g., 
B,V,R,I), cannot be properly adjusted to a model. The fit in- 
cluding the BVRIJHK magnitudes has T e// = 2987 K (M = 
0.12 M Q and is of very low quality (x 2 = 45 for 6 degrees of 
freedom). The fit to BVR alone give T eff = 2880 K, M = 0.11 
M Q and ay 2 = 10.0 for 3 degrees of freedom. In contrast, the best 
fit model obtained from adjusting J, H, and K has ay 1 of 1.92 
for 3 degrees of freedom, an effective temperature of 3225 K and 
a mass of M — 0.172 M©. Note that the mass value is very good 
agreement with the mass estimate derived from the DF00, and 
the effective temperature is surprisingly close to one from the at- 
mospheric transfer models in the previous Section \32] How the 
uncertainty in the stellar mass affects the planet parameters will 
be discussed and accounted for in Section |4] 



4. Star-planet parameters from direct observables 

4. 1 . Using the stellar mass as the input 

Information from the transit light curves together with the 
orbital solution can be combined with the stellar mass (or 
the stellar radius) to fully characterize the bulk properties 
of the star-planet system. These methods have been devel- 
oped over the years by different authors and we suggest read- 
ing ISeager & Mallen-Ornelasl d2003l) . ISouthworthl d2008l) and 
ICarter et aD(l201 ll) for more information on the derivations. This 
subsection describes the relevant relations and the general ap- 
proach used to derive the star-planet parameters using the stellar 
mass as the input (mass-input approach). 



Assuming a circular orbit for the planet, the transit light 
curve alone allows one to obtain a direct measurement of the 
mean stellar density p», c , rc . Given a fully Keplerian solution from 
the Doppler data, the stellar density also depends on the eccen- 
tricity of the orbit and can be written as 



p* =p* 



f VT 



1 + e sin to 



(1) 



Because no prior information is required from the orbital fit, 
p^circ is a quantity typically provided by the studies analyz- 
ing transit light curve s of GJ 1214 (e.g., ICarter et al.l 120111: 
iKundurthv et alJ 1201 lh . These two studies present the analy- 
sis of new light curves and combine the m with previous li ght 
curves from lCharbonneau et all {2009) and lSada et afl d2010h . A 
weighted mean of p*, cl > e from these two studies will be used in 
all that follows. The detailed derivation of Eq.Q]was first given 
bv ISeager & Mallen-Ornelasl (120031) for circular orbits, while 
Carte r et al.l ( 201 ll) included the dependence on the eccentricity. 

With the mean stellar density from Equation[T]and the stellar 
mass from the empirical calibrations discussed in Section 13.31 
one can trivially derive the radius of the star. Using this radius 
and the transit depth R p /R* (again, a direct observable from the 
light curves) one then obtains the radius of the planet. The com- 
bination of the minimum mass (from RV) and the inclination 
(from transit light curve) provides the planet's true mass, which 
is then combined with the planet radius to finally derive the mean 
planet density. 

To obtain realistic a posteriori distributions, one needs to 
assume realistic distributions for the input parameters of the 
model. The complete list of input parameters used at this point 
are given in Table [6] For all the values borrowed from the lit- 
erature, a Gaussian distribution N <x] is assumed where p is 
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Table 6. Input parameters used to compute the Monte Carlo dis- 
tributions of the star-planet parameters 



Parameter 


Distribution 


Expected 


Standard 


Jxcl. 


name 


type 


value 


deviation 




M. [M e ] 


Gaussian 


0.176 


0.009 


1,2 


r -In 
P,,circ [gem 3 ] 


Gaussian 


23.695 


1.7 


3,4 


a/R» 


Gaussian 


14.62 


0.3 


3,4 


R P /R, 


Gaussian 


0.01178 


0.001 


3,4 


Inc. [deg] 


Gaussian 


89.19 


0.5 


3,4 


Period [days] 


Bayesian 


1.580400 


1.4 10 -5 


1 


K[ms-'] 


Bayesian 


10.9 


1.6 


1 


J-o [deg] 


Bayesian 


210.9 


6.2 


1 


e cos w 


Bayesian 


-0.033 


0.055 


1 


e sinw 


Bayesian 


-0.044 


0.090 


1 



Notes . (A) (1) This work (2)lDelfosse et alJd2000t) , (3) lKundurthv et alj 
d201lh . (4) ICarter et al.1 d201 lh 



the preferred value and <x is its published uncertainty. For the 
parameters derived from the orbital solution in Section [23] we 
directly draw samples from the MCMC distributions generated 
during the orbital analysis. The uncertainty in the stellar mass 
has to be accounted for at this point. Given the precision in the 
distance and in the J, H and K photometry, the major source of 
uncertainty in the mass is due to the actual accuracy of the DF00 
calibration. This accuracy is not very well known, but other stud- 
ies suggest that it should be correct at the 5% level, which is the 
relative uncertainty we will use for the stellar mass. 

The process of solving for all the star-planet parameters is 
done for the 10 6 synthetic input sets. The result is a numerical 
representation of the empirical probability distributions for the 
derived star-planet properties. In Figure [7] we show the obtained 
distributions concerning the planet properties only (mass, radius, 
density and surface gravity). The surface gravity of the planet 
can also be de rived from observa bles only using the prescrip- 
tions given by ISoumworthl d2008h (combination of light curve 
parameters and parameters from the orbital solution). Table [7] 
provides the expected values and standard deviations of the dis- 
tributions for each derived star-planet parameter. 

We find that the mean density of GJ 1214b has to be smaller 
than 2.40 g cirT 3 with a 95 % confidence level (c.l.), with 1.69 
g crrT 3 being the expected value of the distribution. An Earth- 
like density (e.g., p > 5.5 g cm 3 ) is thus ruled out at a > 99.9 % 
c.l. As illustrated in Fig. [7] this density range confirms that the 
planet is more si milar to a small ve rsion of Neptune, or a water 
dominated body dBerta et al.ll2012l e.g., ocean planets, ), rather 
than a scaled-up version of the Earth with a solid surface. This 
information by itself does not solve the question about the atmo- 
spheric composition discussed in the Introduction, but, at least, 
eliminates the uncertainty in the distance of the star in modeling 
of the possible planetary structure. 

We want to stress that, thanks to the new distance estimate, 
we now find a remarkably good agreement of the derived stellar 
radius compared to the one obtained from the SED fit in Section 
1121 

As a final test, we reproduced the star-planet parameters us- 
ing the stellar radius instead of the stellar mass. This consists 
of doing the following: the stellar radius (SED fit) combined 
with the mean stellar density (light curve plus RV) gives the 
stellar mass. Then the stellar mass combined with the RV ob- 
servables and the inclination (light curve) provides the planet 
mass. In parallel, the stellar radius with the transit depth (light 
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Planet surface gravity [m s" 2 ] Planet density [g cm" 3 ] 



Fig. 7. Final probability distributions for most significant bulk 
properties of GJ 1214b. The mean densities of a few represen- 
tative Solar system objects are marked as black vertical bars for 
reference. 

Table 7. Suggested new star-planet parameters as derived using 
the stellar mass as the only input parameter. Note that some of 
the distributions in Fig.[7]are strongly asymetric. 



Parameter Maximum Expected Standard 

name probability 1 value deviation 



Star mass [M Q ] 


0.176 " 




0.0087 


Star radius [Re] 


0.213 


0.211 


0.011 


Star surface gravity [m s" 1 ] 


1032 


1109 


220 


logg [g in cm s~']' J 


5.01 


5.04 


0.07 


Star mean density [g cm 4 ] 


26.36 


27.7 


8.9 


Planet radius [R ffi ] 


2.80 


2.72 


0.24 


Planet mass [M e ] 


6.26 


6.19 


0.91 


Planet surface gravity [m s~ 2 ] 


7.66 


7.68 


1.19 


Planet mean density [g cirT 3 ] 


1.56 


1.69 


0.61* 


Planet T C ,[K] (Bond albedo = 0) 




576 


14 


Planet T fc/ [K] (Bond albedo = 0.75) 




407 


11 



Notes. (,) Maximum of the marginalized posterior distribution ( "> Used 
as input {b} From light curve and Doppler analysis (no spectral adjust- 
ment) {c} p p < 2.40 g cm" 3 with a 95% c.l. 



curve) gives the planet radius, which combined with the mass 
finally provides the planet density. This approach, however, has 
a non-obvious drawback. The mass obtained through the stellar 
density formula depends as the third power on the stellar radius 
(M* oc /? 3 ), making the stellar mass determination very sensitive 
to small radius changes. Going the other way around (starting 
from the mass, derive the radius), the stellar radius depends as 
Mr'^ 3 , resulting in a much weaker dependence. For example, a 
5% uncertainty in the mass translates to a 1 .6% uncertainty in the 
radius. On the other hand, a 5% uncertainty in the radius trans- 
lates to a 15% uncertainty in the mass. When the uncertainty in 
the stellar density is also included, this difference is even more 
pronounced. Therefore, we consider the mass input method a 
much more robust way of providing a consistent picture. The 
updated values using the mass input approach for the star-planet 
parameters are in Table Q and the marginalized distributions for 
the relevant planet parameters are depicted in Fig. [7] 
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5. Conclusions 

The metallicity of GJ 1214 (and of M dwarfs in general) and its 
effects on the observable fluxes (e.g., optical fluxes) still require 
a better understanding, both observationally and theoretically. 
Even though spectroscopic methods in the nIR are getting bet- 
ter at determining the metal content of cool stars, precise direct 
measurements of their distances are still required for the charac- 
terization of exoplanets around M dwarfs. At least for the stel- 
lar masses, the model predictions from the near-infrared fluxes 
seem to be in good general agreement with current empirical 
relations derived from measured masses of M dwarfs in bina- 
ries (e.g., DFOO). This is fortunate because the stellar mass is 
the only input parameter required to derive the star-planet bulk 
properties from direct observables. The fit of the stellar spec- 
tral energy distributions to the BT-Settl-2010 model grid also 
provides consistent estimates in the stellar properties if infrared 
photometry is used. Even with the remaining ambiguities in the 
stellar parameters, both approaches (using a stellar mass from 
empirical relations, or deriving a radius from the absolute fluxes) 
lead to consistent results for the star-planet bulk properties, but 
the mass input approach is preferred due to the lower sensitivity 
of the method to uncertainties in the input parameters. Several 
studies have been published in recent years trying to better char- 
acterize the GJ 1214 star-planet system. The major source of the 
reported uncertainties came from directly observable quantities 
that we have improved, collected, and combined here: trigono- 
metric parallax and corresponding absolute fluxes, improved RV 
measurements, inclusion of all of the transit observations in the 
derivation of the orbital solution, additional infrared flux mea- 
surements, and light curve observables derived from photomet- 
ric follow-up programs. 

We now find remarkable agreement of the derived star prop- 
erties obtained by comparing the JHK fluxes to up-to-date at- 
mospheric and evolutionary models. All previous studies had to 
invoke some mechanism (e.g., spot coverage) to justify the mis- 
match between the predicted versus observed properties of the 
star. This alone highlights the importance of obtaining direct and 
accurate distance measurements of low mass stars. 

At this time, the quantity that most requires further improve- 
ment and/or independent determination is the radius of the star. 
To our knowledge, this can be observationally achieved with 
3 different methods: 1) additional RV measurements to fur- 
ther constrain the orbital eccentricity, 2) direct measurement 
of the stellar diamete r using optical/nIR interferometry (e.g., 
Ivon Braun et al.l 1201 2l) . or 3) detection of the secondary tran- 
sit (whose instant strongly depends on e and a>). We note that 
GJ 1214 is faint at optical wavelengths compared to other stars 
observe d by high precision RV instruments (e.g., HARPS or 
HIRES iBonfils et alJ 1201 it IVogt et all l2010h . In the case of 
HARPS, integrations of 45 min were required to obtain a pre- 
cision of ~3.4 m s _1 , so a refinement of the orbital eccentric- 
ity through additional RV measurements is time-consuming but 
quite possible. On-going space-based photometric observations 
in the mid-infrared (e.g., using Warm Spitzer/NASA) should be 
able to detect the secondary transit soon and pin down the or- 
bital eccentricity to greater precision. We hope that the updated 
orbital solution provided here facilitates the task of finding such 
a secondary transit. 
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Table 8. Transit time observations used in the orbital fitting 
BJD Uncertainty Source 



(days) 




(days) 




2455307 


892689 


0.000263 


1 


2455353 


724652 


0.000311 


1 


2455383 


752334 


0.000264 


1 


2454980 


748942 


0.000417 


1 


2454983 


909686 


0.000401 


1 


2454964 


944935 


0.000788 


2 


2454980 


748976 


0.000264 


2 


2454983 


909689 


0.000228 


2 


2454999 


713690 


0.000253 


2 


2454980 


748570 


0.000150 


3 


2454983 


909820 


0.000160 


3 


2454988 


650808 


0.000049 


3 


2455002 


874670 


0.000190 


3 


2455269 


962990 


0.000160 


3 


2455288 


928200 


0.001100 


3 


2455296 


830130 


0.000230 


3 


2455315 


794850 


0.000230 


3 


2455315 


794693 


0.000080 


3 


2455318 


955230 


0.000170 


3 


2455353 


723870 


0.000180 


3 


2455356 


884950 


0.000150 


3 


2455364 


787000 


0.000150 


3 


2455375 


849970 


0.000130 


3 


2455383 


752050 


0.000130 


3 


2455391 


654105 


0.000059 


3 


2455315 


793430 


0.000420 


4 


2455345 


821260 


0.000140 


4 


2455345 


821330 


0.000370 


4 


2455353 


723320 


0.000360 


4 


2455364 


786690 


0.000290 


4 



Notes. (l)(Kundurthv et al 



analiz ed light curves from 
120111) (4) dSadaetalJl2O10h 



.20111) (2)dKundurthy et al]|201l|) using re- 
Charbonn eau et alj d2009l) (3) dCarter et alj 
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